Surface morphology and magnetic anisotropy in (Ga,Mn)As 
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Atomic Force Microscopy and Grazing incidence X-ray diffraction measurements have revealed the presence of ripples 
aligned along the [110] direction on the surface of (Ga,Mn)As layers grown on GaAs(OOl) substrates and buffer layers, 
with periodicity of about 50 nm in all samples that have been studied. These samples show the strong symmetry 
breaking uniaxial magnetic anisotropy normally observed in such materials. We observe a clear correlation between the 
amplitude of the surface ripples and the strength of the uniaxial magnetic anisotropy component suggesting that these 
ripples might be the source of such anisotropy. 

PACS numbers: 75.50.Pp, 75.30.Gw, 07.79.Lh, 61.05.cp 



The ferromagnetic semiconductor (Ga,Mn)As has been the 
subject of intense research due to its potential application in 
spintronic devices and quantum computation 12 . The ferro- 
magnetic interaction between the local Mn moments is medi- 
ated by the itinerant holes giving rise to a strong sensitivity 
of the magnetic properties on the carrier density and lattice 
strain 3 4 . Developing an understanding of how these interac- 
tions give rise to the magnetic properties is essential in order 
to realise the applications of the (Ga,Mn)As compound. The 
p-d Zener model gives a satisfactory explanation for several 
of the properties. For instance, (Ga,Mn)As grown epitaxi- 
ally on GaAs experiences a compressive which gives rise to 
a magnetic anisotropy with the easy axis in the plane of the 
layer, consistent with the predictions of the model 34 . In ad- 
dition, the in plane magnetic anisotropy is found ubiquitously 
to consist of a cubic contribution, favouring easy axes along 
the [100] directions, and a uniaxial component favouring one 
of the [110] directions 56 . The presence of the cubic com- 
ponent is consistent with the symmetry of the GaAs lattice, 
but the microscopic origin of the uniaxial component has not 
yet been explained. Several theoretical studies have invoked 
the concept of a uniform anisotropic strain o f ord er 10 ~ 4 to 
account for the uniaxial anisotropy component 4 5 7 , but no di- 
rect experimental evidence of such a lattice distortion has been 
reportecP. 

In this paper we show experimental evidence, from Atomic 
Force Microscopy (AFM) and Grazing incidence X-ray 
diffraction (GID) measurements, of the presence of periodic 
ripples in (Ga,Mn)As layers grown on GaAs. The analysis 
of AFM and GID measurements have evidenced a period of 
these ripples of about 50 nm. We observe a clear correlation 
between the magnitude of the uniaxial anisotropy and the am- 
plitude of the surface ripples. 

(Ga,Mn)As films with thicknesses (t), 5nm, Inm and 
25nra, and Mn concentrations, 6% and 12%, have been 
deposited on GaAs (100) substrates by low temperature 
(~200°C) molecular beam epitaxy on lOOnm thick low tem- 
perature or high temperature (~580°C) GaAs and AlGaAs 



buffer layers. Atomic force microscopy (AFM) images were 
obtained by using an Asylum Research MFP-3D atomic force 
microscope. Before performing the measurements the sam- 
ples were cleaned with a 1:3 HCkF^O solution to remove the 
superficial oxide layer. Measurements on uncleaned samples 
have shown the same surface structures but with additional 
surface contamination indicating that the cleaning was not 
modifying the surface morphology. The images were taken 
on lxl fim 2 areas along two different crystallographic direc- 
tions, [110] and [110]. Grazing incidence X-ray diffraction 
(GID) measurements were performed at the beamline ID01 of 
the ESRF in Grenoble (France) using a photon energy of 6.5 
keV and a linear detector mounted perpendicular to the sample 
surface. Magnetometry measurements on the as-grown films 
were carried out using a Quantum Design Superconducting 
Quantum Interference Device (SQUID) magnetometer. 

Fig.[T]shows the remnant magnetisation measured along the 
[110] and [110] directions for one of our samples. Consis- 
tent with previous studies of unannealed (Ga,Mn)As layers 
we find that all samples within our study have a dominant 
in plane [110] uniaxial magnetic easy axis at high tempera- 
ture, with competition between in plane [110] uniaxial and 
in-plane [100] bi-axial easy axes at low temperatures. We 
consider a phenomenological description of the free energy in 
terms of cubic (Kc) and uniaxial (Ku) anisotropy constants: 
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FIG. 1 . Temperature dependence of the remnant magnetisation along 
[110] and [110] direction. 
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FIG. 2. AFM image of the sample Mn 381 a) along [110] direction, 
andb) along [110] direction. 

E = K c sin 2 (0) cos 2 (0) + K v sin 2 (0) - MsHcos(9 - O ), 
where 6 and 0q measure the direction of the magnetisation 
and the external magnetic field with respect to the [110] crys- 
tal direction. By minimising this expression for the energy it 
can be shown that the relative strength of the uniaxial to cu- 
bic anisotropy Ku /Kc = cos 2[tan _1 (M 1 Y /Mno)], where 
M x y and Mn represent the components of the magnetisa- 
tion measured along [110] and [110] respectively. We find that 
the temperature-dependence of the remnant magnetization af- 
ter field cooling in all the studied samples is generally well- 
described by this phenomenological model. 

Fig. [2] shows AFM images along the [110] and [110] direc- 
tions for a 7nm thick (Ga,Mn)As sample containing 6% Mn 
(sample Mn381). The images clearly reveal the presence of 
ripples aligned along the [110] direction. Ripple patterns such 
as these are known to form during the MBE growth of low 
temperature (< 260° C) GaAs in the presence of excess As. 
The mechanism causing the ripples to form is the suppres- 
sion of diffusion of adatoms at step edges. The excess As in 
ref P was believed to act as a surfactant aiding the diffusion 
of adatoms along the step edges. For our samples, which are 



sample t , Mn% 


buffer 


T 


RMS 


K u /K c 


Mn 352 


5,6% 


AHT 


49.7 ± 12.7 


0.35±0.07 


0.241 


Mn536 


5,6% 


BHT + BLT 45.5 ± 11.3 


0.31±0.04 


0.157 


Mn 396 


7, 6% 


BHT 


47.7 ± 8.2 


0.48±0.06 


0.498 


Mn381 


7,6% 


BHT 


41.9 ±6.4 


0.38±0.04 


0.325 


Mn 555 


7, 6% 


AHT 


47.7 ± 6.3 


0.3±0.1 


0.158 


Mn535 


7, 6% 


BLT 


45.4 ± 10.3 


0.34±0.04 


0.089 


Mn 394 


7, 6% 


AHT 


59.8 ± 13.5 


0.35±0.05 


0.249 


Mn 437 25, 12% 


ALT 


51.5 ±8.6 


0.36 ± 0.06 


0.113 


Mn 438 25, 12% 


ALT 


61.7 ± 11.4 


0.37 ±0.07 


0.149 


Mn 554 


25, 6% 


BLT 


56.0 ± 16.9 0.35 ± 0.04 


0.039 


Mn 499 


25, 6% 


AHT 


47.8 ± 11.8 


0.31±0.05 


0.018 


Mn 330 


25, 6% 


AHT 


65.9 ± 11.3 0.31 ±0.09 


0.107 


Mn 490 


25, 6% 


AHT 


52.5 ± 8.3 


0.28 ± 0.03 


-0.005 



TABLE I. For each sample the thickness (t) in nm, the Mn concen- 
tration (Mn %), the buffer layer (A: GaAs; B: AlGaAs; HT: high 
temperature; LT: low temperature), the period (r) in nm, the RMS 
roughness in nm and the magnetic anisotropy coefficient {Ku / Kc), 
are shown. 



not grown with a significant excess of As, we believe that Mn 
deposited at the surface may be playing a similar role as a sur- 
factant. Similar ripple patterns are observed in the images ob- 
tained for all the samples in our study. We have used Fourier 
analysis to extract quantitative information from these images. 
We can define the power spectral density as p s = | (/> s | 2 , where 
(j) s is the Fourier amplitude corresponding to a (dimension- 
less) frequency f s = s — 1. 

The mean frequency, f(y), per line is then obtained as the 
average of the frequencies of all Fourier harmonics (along a 
given scan line) weighted with the corresponding power spec- 
tral densities: /(2/) = [T,%l( s - l )Ps]/[YTs=iPs\. Then, 
the global 'effective' mean-frequency (/), characterising the 
oscillations in the whole image, is given by the mean f(y) of 
each line along all horizontal scan lines and the associated er- 
ror is its standard deviation. From the / and L (the width of 
the image) we obtain the 'effective' period of the oscillations: 
r = L/f. Such an effective period is used to give a quantita- 
tive measure of the periodicity of the ripples observed on the 
AFM image. In addition, the root mean square (RMS) rough- 
ness of each image in the direction perpendicular to the ripples 
is calculated by taking the standard deviation of the height for 
each line scan and averaging this over the whole image. 

The ripples are found to be present in all of the (Ga,Mn)As 
samples studied. In Fig.[3]we plot Ku / Kc from SQUID mag- 
netometry measurements at 2K versus the RMS roughness, 
and in Table [I] we tabulate the data for each sample. In order 
to distinguish the effect of the roughness from effects which 
may alter the absolute value of Kc we compare the samples 
within groups of similar thickness, nominal Mn concentration 
and Curie temperature. We find that for each group of sam- 
ples there is a clear relation between the ratio of the magnetic 
anisotropy constants Ku/Kc and the RMS roughness, in par- 
ticular the influence of the uniaxial anisotropy increases as the 
RMS roughness increases. 

In order to obtain information on the periodicity of the 
nanostructures over the whole sample surface, grazing inci- 
dence X-ray diffraction (GID) measurements were also per- 
formed on one sample (Mn 381). In the GID geometry, both 
the angle of incidence {on) and the exit angle (a/) of the x- 
ray beams are close to the critical angle(a c ) of total external 
reflection (0.38° for the energy used) thus leading to a large 
footprint illuminating the sample surface, after the aperture 
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FIG. 3. The magnetic anisotropy coefficient, Ku/Kc as a function 
of the total RMS roughness for samples with 6% of Mn. 
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FIG. 4. Transverse GID scans for sample Mn381 collected along 
two perpendicular directions at m = 0.35°. While the scans with 
the primary beam directed across the ripples (a) show two symmetric 
satellites denoting a periodicity of 48 nm, such satellites are absent 
in the scans along the ripples (b). (c) Profile from data in panel (a) 
and a fit assuming Gaussian-shaped ripples are shown by circles and 
solid line, respectively, (d) Sketch of the dimensions of the ripples. 

of horizontal slits, up to a few millimeters. The difference in 
the azimuthal directions of the beams was set close to 20 b, 
where 0b is the Bragg angle 10 . In particular, two in-plane 
diffraction vectors h = 220 (perpendicular to the ripples), and 
h = 220 (parallel to the ripples) were investigated and, for 
each diffraction, we measured two reciprocal- space distribu- 
tions of the scattered intensity in the reciprocal planes Q a Q z 
and QaQz, where Q a Q z is the component of the scattering 
vector Q perpendicular to the surface, and Q a is the angu- 
lar component of Q perpendicular to h. The angular map in 
Fig. [4] (a) taken in the diffraction 220 exhibits two side max- 
ima (intensity satellites). Since Qa in this map is parallel to 
[110], i.e. perpendicular to the ripples, these maxima stem 
from the periodicity of the ripples and the satellite distance 
is inversely proportional to the ripple period r. In turn, the 
width of the satellites depends on the root-mean- square de- 
viation of the ripple spacing a T . From the intensity map we 
extracted the intensity integrated over Q z and compared it to 
a theoretical model. In the model, we assumed that the ripple 
sequence obeys the standard one-dimensional short-range or- 
der modeP^ and that the shape of a cross-section of a single 
ripple is Gaussian: h{x) = hoex.p(—x 2 /R 2 ). 

From Fig.|4jc) comparing the measured integrated intensity 
to the simulations it follows that the model well describes the 
actual ripple structure in excellent agreement with the AFM 
data. From the fit we determined the mean ripple period 
r = (48±5)nm, a T = (14±3)nm, and the shape parameters 
R = (19 ± 5)nm and ho = (0.25 ± 5)nm. In contrast, the 
angular map in the other diffraction in Fig.|4jb) does not show 
any satellites. In this map, the Q a axis is parallel to the [110] 
direction, i.e. along the ripples. The surface morphology does 
not exhibit any pronounced space frequency in this direction 



so that no side maxima is observed. From the fact that the in- 
tensity distribution along the Q a axis in Fig. [4] (a) is perfectly 
symmetric we can argue that the component of the strain par- 
allel to the ripples is either zero or it does not depend on the 
real-space coordinate along [110], i.e. perpendicular to the 
ripples. A more substantial GID-based analysis on a variety 
of (Ga,Mn)As samples will be needed in order to possibly re- 
late the formation of the periodic structures with the presence 
of uniaxial strain, that may lead to a deeper understanding of 
the origin of the magnetic anisotropy; this will be the subject 
of further work. 

In conclusion, AFM and X-ray GID measurements have re- 
vealed the presence of ripples aligned along the [110] direc- 
tion on the surface of (Ga,Mn)As layers grown on GaAs(OOl) 
substrates and buffer layers. We find that the strength of the 
uniaxial anisotropy of these layers increases with the ampli- 
tude of the ripples. This suggests that the symmetry breaking 
uniaxial anisotropy observed in (Ga,Mn)As might be related 
to these symmetry breaking structural features. 
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